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Abstract 


The customary perspective to reason about epidemic mitigation in temporal networks 
hinges on the identification of nodes with specific features or network roles. The ensuing 
individual-based control strategies, however, are difficult to carry out in practice and 
ignore important correlations between topological and temporal patterns. Here we 
adopt a mesoscopic perspective and present a principled framework to identify collective 
featnres at multiple scales and rank their importance for epidemic spread. We use 
tensor decomposition techniqnes to build an additive representation of a temporal 
network in terms of mesostructures, snch as cohesive clusters and temporally-localized 
mixing patterns. This representation allows to determine the impact of individual 
mesostructures on epidemic spread and to assess the effect of targeted interventions 
that remove chosen structnres. We illustrate this approach using high-resolution social 
network data on face-to-face interactions in a school and show that our method affords 
the design of effective mesoscale interventions. 


Keywords: complex networks, temporal networks, tensor decomposition, epidemic spread, 
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Many natural and artificial systems are ade¬ 
quately represented in terms of interaction net¬ 
works between their components^. The network 
representations has led to key insights on the re¬ 
lation between the structure of these systems and 
the dynamics of diverse processes that take place 
over the network, such as information diffusion, epi¬ 
demic spread, and much morel^. In recent years, 
the availability of time-resolved network data has 
pushed network science beyond the static graph 
representation and has prompted new research on 
understanding and modeling time-varying graphs, 
commonly referred to as temporal networkJ^. This 
has spurred an intense activity on studying how 
dynamical processes are affected by the temporal 
evolution of the network on which they unfold. 


focusing on heterogeneous distributions of inter¬ 
event times (“burstiness”), heterogeneous activity 
distributions, causality constraints, etc.^^^^i^. High- 
resolution social network data^^^^^, in particular, 
have opened the way to richer individual-based 
models of epidemic spread^i^l^IHISl and have raised 
new questions on designing control strategies for 
epidemic processes in temporal networkJ^^HUM] 
However, micro-interventions that operate on indi¬ 
vidual nodes, such as targeted vaccination, are dif¬ 
ficult to implement, both for lack of high-resolution 
data in general cases and because action at the 
level of individuals is subjected to many infor¬ 
mational and operational constraints. Therefore, 
it is natural to consider interventions that target 
important collective patterns and structural units 
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of temporal networks. Since network communi¬ 
ties, i.e., cohesive clusters, have long been recog¬ 
nized as key structural units in the architecture 
of complex network^^, it is natura l to t hink of 
interventions at the community levell^^l^Sl Time- 
varying networks, however, can exhibit a richer 
range of activity/connectivity patterns than static 
networks. Defining and detecting temporal net¬ 
work structures at the intermediate (“meso”) scale, 
as well as understanding their relevance for dy¬ 
namical processes occurring on the networks, re¬ 
mains a largely open question with important ap¬ 
plications. Early work on the mesoscale structure 
of complex networks has mostly focused on static 
networks^ and on community detection in tem¬ 
poral networksl^. Recently, latent factor analysis 
techniques have been used to simultaneously detect 
topological and temporal activity patterns of time- 
varying network^^. 

Here we tackle the problem of designing interven¬ 
tions that selectively target the mesoscopic struc¬ 
ture of a temporal network, with the goal of control¬ 
ling a dynamical process such as epidemic spread. 
This requires both detecting mesostructures and 
ranking them by their importance for a specific 
dynamical process. To detect mesostructures we 
build on the work of Ref. |3I] and use tensor de¬ 
composition technique^^ to represent a temporal 
network as an additive superposition of mesostruc¬ 
tures, each of which can encode complex structural 
and temporal correlations. To assess the impor¬ 
tance of individual mesostructures, we take advan¬ 
tage of the additivity of the tensor decomposition, 
which allows us to dissect the original network and 
reassemble altered temporal networks where chosen 
mesostructures are selectively removed. By com¬ 
paring the dynamics of a given process on the orig¬ 
inal network and on the altered ones we can then 
determine the specific role of individual mesostruc¬ 
tures and we can rank the effectiveness of interven¬ 
tions that target those structures for removal. 

We illustrate our approach using high-resolution 
social network data on face-to-face interactions in a 
school setting and investigating the importance of 
individual mesostructures for simple epidemic pro¬ 
cesses. Schools are actually an interesting context 
for models of epidemic spread, as they are thought 
to play an important role in the community spread 
of infectious disease^^^HlZI Macro-scale interven¬ 
tions such as school closure are considered a vi¬ 


able strategy for epidemic miti gation !^ but come 
with steep socio-economic cost^^^Eo!^ thus calling 
for more targeted approaches. We show that our 
approach identifies effective strategies for epidemic 
mitigation that involve the targeted removal of in¬ 
dividual mesostructures. Remarkably, we find that 
the most important mesostructures for epidemic 
spread correspond neither to the strongest struc¬ 
tures nor to cohesive structures that can be easily 
detected by community detection methods. 

Mesoscale structure of temporal networks 
Latent factor analysis of temporal networks 

We first summarize a general framework that we in¬ 
troduced in ref. ED to expose the mesoscale struc¬ 
ture of a temporal network. The approach com¬ 
bines a tensor representation of a temporal net- 
worklSI and a dimensionality reduction technique 
based on tensor decomposition!^. 

We consider an undirected and unweighted tem¬ 
poral network, hence the state of the network at 
time t can be represented by a binary-valued ad¬ 
jacency matrix M(t) € where the matrix 

entry (t) indicates the status of link i-j at time 
t, and N is the number of nodes of the network. 
The matrices describing the state of the network at 
different times can be combined into a 3-mode ten¬ 
sor, T G ^nxNxl^ where the first two dimensions 
are the customary node indices of the adjacency 
matrix and the third dimension is a temporal in¬ 
dex. L is the number of network snapshots, each 
for a different point in time. To detect structures, 
i.e., correlated activity/connectivity patterns of the 
temporal network, we set up a latent factor anal¬ 
ysis based on non-negative tensor factorization!^ 
The central idea is to approximate the tensor T 
with a sum T of rank-1 non-negative tensors (see 
Methods): 

R 

( 1 ) 

r—1 

with R < minjfV, S'}. This decomposition al¬ 
lows to describe the original temporal network as 
a purely additive superposition of a chosen num¬ 
ber of activity/connectivity patterns, encoded by 
the individual tensor components. Each compo¬ 
nent Sr is regarded as a mesoscale structure of the 
temporal network and can encode cohesive clus- 
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ters (i.e., communities), temporally-localized mix¬ 
ing patterns, and more. Such a parts-based repre¬ 
sentation is a general and powerful feature of tensor 
decompositions with non-negativity constraints^ 
and greatly enhances the interpretability of individ¬ 
ual components. A similar approach based on non¬ 
negative matrix factorization has been proposed 
to detect overlapping communities in static net- 
workS^. 

The additive decomposition of Eq.[^also suggests 
a natural way to investigate the role played by in¬ 
dividual mesoscale structures, both on the overall 
network architecture and on the dynamics of pro¬ 
cesses occurring over the network, e.g., epidemic 
spread. In particular, we can assess the impor¬ 
tance of individual mesostructures for dynamical 
processes unfolding over the temporal network by 
selectively removing those structures: We can start 
from the decomposition of Eq. remove a cho¬ 
sen term (i.e., a chosen structure), and sum all the 
other terms, obtaining a modified temporal net¬ 
work where the chosen structure has been selec¬ 
tively erased. We can then simulate a given dy¬ 
namical process over the original temporal network 
and over the modified temporal network, and com¬ 
pare the dynamics or the outcome of the process, 
hence learning about the impact of the removed 
structure. 

Case study 

As a case study we consider an empirical tempo¬ 
ral network of human face-to-face interactions in a 
primary school!^, measured by using wearable sen- 
sors^. The dataset describes the interactions (or 
“contacts”) between 231 children aged 6 to 12 and 
10 teachers, organized in 10 classes. It comprises 
the contacts that occurred at the school premises 
during two consecutive days in October 2009, from 
8:30am to 5:15pm of each day. We ch ose this 
dataset because it is publicly available^^^^^ and pro¬ 
vides not only a high-resolution temporal network, 
but also ground truth information on communities, 
namely, class attributes for all nodes. Moreover, 
the school schedule is known and can be used to 
understand the activity patterns that mix multi¬ 
ple classes, as done in Ref. [3TJ In the following this 
will allow us to relate the effect of removing an indi¬ 
vidual component to specific knowledge about the 
social behavior associated with that component. 


The school data are naturally represented as a 
temporal network, where nodes are individuals and 
links are interactions between pairs of nodes. The 
networks snapshots of our dataset are recorded ev¬ 
ery 20 seconds, however the school activity is sched¬ 
uled on a coarser temporal scale. Hence, we aggre¬ 
gate the data over consecutive time intervals of 15 
minutes to obtain unweighted network snapshots: 
for every 15-minute interval, a link is drawn be¬ 
tween two nodes if the corresponding individuals 
had at least one contact during that interval. This 
yields a temporal network represented by binary¬ 
valued valued tensor T with N = 241 and L = 70. 


Following ref. 1311 we carry out the non-negative 
tensor factorization of Eq. and approximate the 
tensor T with the sum T of i? = 14 compo¬ 
nents (see Methods and Supplementary Informa¬ 
tion) : the results of the decomposition are reported 
in Fig. [2 Components through Sg and compo¬ 
nent iSii correspond to the 10 school classes of the 
school. They can be validated using the know class 
attributes for students and teachers, and they can 
also be detected to some extent using community 
detection algorithms on the temporally-aggregated 
network!^. The components 5io, 5i2, 5i3 and S14, 
however, mix individuals from different classes and 
- contrary to the above class components - they 
are only active during lunch breaks. Such compo¬ 
nents cannot be found by community detection of 
the temporally-aggregated networks as they corre¬ 
spond to temporally-localized mixing patterns of 
the temporal network, due to mixing of students 
at lunch time. Tensor decomposition, however, can 
naturally detect these patterns as well as the cus¬ 
tomary cohesive network communities. 


We remark that whereas the original temporal 
network T is unweighted, the approximating net¬ 
work T is in general weighted (see Methods). For 
any given network snapshot, corresponding to a 15- 
minute time interval, we will interpret the weight 
of the edge between two nodes as the cumulated 
time that those nodes spent in contact over that 
15-minute interval. In Table we report the frac¬ 
tion of total tensor weight corresponding to each of 
the i? = 14 components, ||5r|| / ||T||. 
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NODES 


DAY 1 ACTIVITY 


DAY 2 ACTIVITY 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 8:00 10:00 12:00 14:00 16:00 8:00 10:00 12:00 14:00 16:00 


Figure 1: Components obtained by tensor factorization of the school temporal network. 

Each row is associated with one component 5^. = ar o ar o c^ of the tensor decomposition T in Eq. 
(see Methods), with R = 14. Components are ordered top to bottom according to their Euclidean norm. 
Left-hand color-coded matrix: Node weights for each component. Each row corresponds to a vector 
ar of the decomposition, normalized to unity in Euclidean norm for clarity of visualization. Nodes are 
ordered to expose the block structure of the matrix due to school classes. Right-hand plots: Component 
activity as a function of time over the 2-day span of the dataset. Each row correspond to a vector Cr, 
normalized to unity in Euclidean norm. Components Si to Sg and component iSn are school classes, 
while components 5io, iSi 2 , iSia and iSi 4 mix individuals from different classes and are active during 
lunch breaks. 


Mesoscale targeted interventions for epi¬ 
demic spreading 

Intervention design and evaluation 

The decomposition of a temporal network into 
a superposition of components, interpreted as 
mesostructures, makes it possible to pinpoint and 
assess the contribution of individual mesoscopic 
features to the overall structural and functional 
properties of the temporal network. Indeed, the ad¬ 
ditive decomposition of Eq. allows to selectively 
remove a chosen component s by excluding the cor¬ 
responding term, yielding an altered temporal net¬ 
work, 

= E ‘5- (2) 

r^s 

which can be compared to the original network to 
elucidate the specific role of component s. When 
component s can be interpreted in terms of a spe¬ 
cific behavioral pattern (e.g., the lunch breaks of 


Fig. 0 its removal can be regarded as the effect 
of an intervention strategy that selectively targets 
that behavior (e.g., removing lunch breaks from the 
school schedule). 

As a first step in exploring mesoscale interven¬ 
tions, here we focus on the case study of the pre¬ 
vious section and investigate the effect of remov¬ 
ing individual mesostructures on the dynamics of 
simple epidemic processes unfolding over the school 
temporal network. We study how the timing and 
size of the epidemic are affected by the removal of 
individual network components, that is, by differ¬ 
ent targeted interventions aimed at removing them. 

We start with the simple case of a susceptible- 
infected (SI) process. Each network node can be 
in either of two states: susceptible (S) or infected 
(I). A susceptible node in contact with an infected 
one becomes infected with a hxed probability A per 
unit time. Once infected, a node stays indefinitely 
in that state. The system is initialized with all 


4 















































r 

1 

2 

3 

4 

5 

6 

7 

weight fraction 

11.3% 

8.6% 

8.8% 

7.1% 

8.3% 

6.9% 

5.7% 

r 

8 

9 

10 

11 

12 

13 

14 

weight fraction 

6.8% 

7.5% 

7.5% 

4.8% 

6.8% 

5.7% 

6% 


Table 1: Component weight. Fraction of total tensor weight corresponding to each term Sr of the 
decomposition T. 


nodes in the S state, except for a single infected 
node (seed node). The timing and duration of con¬ 
tacts between nodes is described by the temporal 
network T, which is a sequence of weighted net¬ 
works: for each network snapshot, the weight of an 
edge between two nodes gives the total duration of 
the contacts between those nodes during the cor¬ 
responding time interval. Although the SI process 
does not describe any realistic infectious disease, it 
represents a paradigmatic dynamical process and 
it is frequently used to investigate the structural 
properties of temporal networks. 

To investigate the effect on epidemic spread of 
removing a chosen mesostructure Sr, we simulate 
the epidemic process on the altered network 
and compare its dynamics to that observed on the 
full temporal network T. We quantify the effect of 
removing component r by measuring the epidemic 
delay ratio® 



where tj is the half-infection time on T when the 
seed is node j, and is the same quantity for the 
altered temporal network The half-infection 
time is measured from the first time the seed node 
infects another node, and the average (•) is com¬ 
puted over all possible seed nodes j and different 
starting times for the SI process (see Supplemen¬ 
tary Information for details). 

We also study the impact of mesoscale fea¬ 
tures on the dynamics of a more realistic epidemic 
process, the Susceptible-Infected-Recovered (SIR) 
model. In this model, S ^ I transitions occur as 
in the SI case when a susceptible node is in contact 
with an infected one, with probability A per unit 
time. Infected nodes recover with a constant rate 
P-, and Recovered (R) nodes no longer take part in 
the epidemic process. To assess the effect of remov¬ 
ing a given mesostructure r we compare the size of 


the epidemics at the end of the process, measured 
as the final number of nodes that are in states I 
or R. We compute the ratio p between the average 
size fT" of the epidemic on the altered network 
and the average size If of the epidemic on the full 
network 'T : 

= 1)1. (4) 

The average is computed over all possible seed 
nodes and over different seeding times near the be¬ 
ginning of the data set (see SI for details). Note 
that, due to the hnite temporal span of our net¬ 
work data, the epidemic process may not be over 
by the end of the simulation (i.e., some nodes may 
still be in the I state). In this case Eq. only 
yields an estimate of how much the epidemic has 
been mitigated during the time interval covered by 
the data. 

Targeted interventions: case study 

We start from the decomposition of the school tem¬ 
poral network described in Fig. and study the ef¬ 
fect of removing individual mesostructures on the 
dynamics of an SI process. We consider several 
values of the transmission rate A: here we present 
results for A = 0.7 and in the Supplementary Infor¬ 
mation we show that qualitatively similar results 
hold for other values of the parameter. Figure 
shows the delay ratio Tr of Eq. obtained by se¬ 
lectively removing each of the components Sr, with 
r = 1,... ,14: that is, comparing the dynamics 
of the SI process over each of the altered networks 
with the dynamics observed for the full tem¬ 
poral network T. Since the altered networks have 
by dehnition a smaller total weight ||T’’|| than the 
full network, we expect epidemic spread to be mit¬ 
igated for the altered networks. To assess whether 
a given component r plays a structural role for epi¬ 
demic spread that goes beyond its weight, it is thus 
important to compare the delay ratio Tr observed 
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on removing that component with a suitable null 
model. Such a null model is obtained by measur¬ 
ing the delay ratio for several stochastic realizations 
of altered temporal networks obtained from T by 
removing, at random, a fraction of weight equal to 
that of component r (as given by Table [^. The 
distributions of Tr for the null models are shown by 
the boxplots of Fig. 

0.5 

0.4 

0.3 

Tr 

0.2 
0.1 

^■^ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 
index r of the removed component 

Figure 2: Epidemic delay ratio as a function 
of removed component. The epidemic delay ra¬ 
tio Tr of Eq.j^is plotted (solid circles) as a function 
of the removed network component and compared 
to the null model described in the main text (box- 
plots). In each boxplot, the central band is the 
median, the bottom and top of the box are the 
first and third quartiles. Whiskers extent to the 
most extreme data point within 1.5 IQR of the in¬ 
ner quartiles. Points beyond whiskers are outliers. 



We observe that the delay ratio obtained by se¬ 
lectively removing a component (solid circles) is al¬ 
most always larger than the one obtained by re¬ 
moving links at random (boxplots): Such targeted 
interventions have a stronger effect than random 
removals of links. For the first 9 structures, the 
delay ratio is positively correlated with the com¬ 
ponent weight of Table as we could naively ex¬ 
pect. Strikingly, the removal of components 12 
or 14 slows down the epidemic spread consider¬ 
ably. These two mesostructures, despite carrying 
a smaller total weight than most other components 
(see Table , thus appear to play a more impor¬ 
tant role for SI spreading. As shown in Fig. [^both 
components 12 and 14 mix nodes from different 


classes and have activity concentrated during the 
lunch break of the first day. The other components 
correspond either to individual classes (components 
1 through 9 and 11) or to components, such as 10 
and 13, with little activity on the first day. Overall, 
the SI process is slowed down the most by removing 
comparatively weak structures (in terms of weight) 
that mix classes and occur early in the data. We 
remark that the most effective intervention for miti¬ 
gating epidemic spread thus involves the removal of 
seemingly minor mesostructures that correspond to 
complex correlations between link activity and net¬ 
work structure, and cannot be uncovered by stan¬ 
dard community detection approaches. 

We now turn to the case of an SIR process and 
investigate the effect of mesostructure removal on 
the dynamics of an SIR epidemic over the tempo¬ 
ral network. We compute the epidemic size ra¬ 
tio px,tj.ir) of Eq. i for a wide range of values of 
the parameters A and p, and for all components r. 
The results are shown in Fig for selected com¬ 
ponents and in the Supplementary Information for 
all the other components. We notice that for large 
A values and small /i values (bottom right of the 
heat maps) the epidemic size ratio ~ 1: 

in this region the epidemics spreads and finishes 
fast, hence mesostructure removal cannot mitigate 
its size. The parameter region where mesostruc¬ 
ture removal can have a strong effect, therefore, 
is the arc-shaped region visible, e.g., in the heat 
map for r = 1 (lighter blue/green). The removal 
of class components r = 1 and r = 2 mildly miti¬ 
gates the epidemic for a broad range of parameter 
within that region (almost no effect is obtained for 
the other class components, as shown in Fig. S4 
of the Supplementary Information). Removing the 
class component I has a comparatively higher ef¬ 
fect because of the higher overall weight of that 
component (Table [^. 

On taking a closer look at Fig. we notice that, 
depending on the removed component, significant 
epidemic mitigation can be achieved in two main 
parameter regions within the arc. The first region 
falls within the black contour (high values of A and 
p) and corresponds to parameter values yielding 
epidemics that finish by the end of the second day. 
In this region the epidemic can be strongly mit¬ 
igated (pA,/i — 0.5) by removing mesostructures 
that involve mixing of classes on the first school 
day, namely r = 12 and r = 14. This is con- 
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Figure 3: Epidemic size ratio as a function of SIR model parameters. Each heat map corre¬ 
sponds to a targeted intervention that selectively removes one component. For each removed component 
r the heat map shows the epidemic size ratio p\^^{r) as a function of the SIR parameters A and p. 
P\,^i{f) = 1 indicates that the intervention does not affect epidemic size. The white area is the region 
where the epidemic dies out, i.e., it fails to affect more than 1% of the network nodes. The region inside 
the black contour line corresponds to parameter values such that the SIR epidemic finishes within the 
finite span of the school dataset (2 days), both for the full and the altered temporal networks. That 
is, for those parameter values the epidemic size ratio is not affected by the finite temporal span of the 
dataset. 


sistent with the results of Fig. for the SI case. 
Conversely, removing the class-mixing components 
r = 10 and r = 13, which are mostly active on the 
second day (Fig. has a negligible effect in the 
same parameter region. For a sample parameter 
choice within this region, in Fig. we illustrate the 
temporal evolution of the epidemics and the effect 
of removing individual components: the epidemic 
on the unmodified network peaks on the first day 
and, consistently with the above remarks, removing 
r = 12 or r = 14 strongly mitigates the fraction of 
affected nodes, whereas the removal of r = 10 and 
r = 13 has a negligible effect. We remark that re¬ 
moving the early class-mixing mesostructures also 
mitigates the SIR epidemic for a wide range of pa¬ 
rameter values outside the black contour, as visible 
in Fig. and that removing these mesostructures 
is more effective, in general, than removing class 
components, even though the latter account for a 
larger network weight. The second parameter re¬ 


gion of interest lies in the lower part of the arc (low 
values of A and p) and is visible in Fig. [^for r = 10 
and r = 13. This region corresponds to parameter 
values for which the SIR process is slower and the 
epidemic does not finish over the 2-day temporal 
span of the data. In this case, removing either of 
the class-mixing mesostructures r = 10,12,13,14 
as a similar and limited effect on the epidemic size 
reached at the end of the second day. 

Discussion 

We have put forward a methodology to system¬ 
atically assess the relevance of mesostructures for 
dynamical processes such as epidemic spread in 
temporal networks. The proposed method only 
uses time-resolved connectivity information with¬ 
out metadata and allows to uncover complex pat¬ 
terns of correlated link activity whose removal has a 
strong impact on simulated epidemic spreading for 
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day 1 day 2 



day 1 day 2 

Figure 4: Fraction of infected nodes as a 
function of time for a sample parameter 
choice. Temporal evolution of the average frac¬ 
tion of infected nodes (top) and of the cumulative 
fraction of infected nodes (bottom) as a function 
of time for a sample parameter choice of A = 0.11 
and p. = 0.16. The curve marked as none corre¬ 
sponds to the spread on the unmodified temporal 
network T. The other curves show the effect of se¬ 
lectively removing components r = 1,10,11,12,14. 
For the chosen parameter values, epidemic spread 
occurs on the first day, hence removing the early 
class-mixing components r = 12 or r = 14 has a 
strong mitigation effect. 


a broad range of model parameters. The removal 
of an individual mesostructure can be regarded as 
a targeted intervention at the mesoscale, such as 
the removal of a specific group of nodes or the sup¬ 
pression of a specific activity pattern (e.g., a lunch 
break). These targeted interventions have a po¬ 
tentially high relevance, as they do not involve the 
system as a whole, nor they need actions at micro¬ 
scopic level of individual nodes. Hence, they might 
strike a balance between cost and effectiveness that 
might make them actionable in a variety of appli¬ 
cations. 


As an illustration of the method, we have consid¬ 
ered a case study where simple but paradigmatic 
spreading processes occur on top of an empirical 
temporal network describing face-to-face interac¬ 
tions in a school. In this case, the most impor¬ 
tant mesostructures for epidemic spread are not 
those involving most of the link activity, but instead 
consist of weaker, temporally-localized mixing pat¬ 
terns corresponding to scheduled social activities. 
It is important to remark that these mesostructures 
cannot by found by means of community detection 
techniques on static networks, and that our meth¬ 
ods reveals them in an unsupervised fashion. We 
also notice that the actionability of an intervention 
that suppresses a given mesostructure hinges on the 
recognizability of that mesostructure in terms of 
relevant contextual information, such as a known or 
understandable grouping of nodes, temporal activ¬ 
ity localized at specific times of a known schedule, 
or a specific relation to other metadata for the con¬ 
text at hand. For example, the selective removal of 
component iSi 2 of Fig. only becomes an action¬ 
able targeted intervention when its temporal activ¬ 
ity profile is identified as a mixing pattern associ¬ 
ated with the lunch break on the first school day. 


Our work is a first step towards the systematic 
design of mesoscale interventions in temporal net¬ 
works. Natural future steps include the use of our 
methodology on different temporal networks, both 
empirical and synthetic, the investigation of the im¬ 
pact of incomplete or noisy data, and the study of 
other dynamical processes. 
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Methods 

A 3-inode tensor T, with entries Tijk, can be ap¬ 
proximated by a sum T oi R rank-1 tensor^^HlIIlll^ 

R R 

T = ^ ^ ar o br o Cr , (5) 

r—1 r—1 

where each component tensor Sr is the outer prod¬ 

ucts of three vectors, Sr = o br o Cp. The 
Frobenius norm of the difference between T and 
T, e = ||T — T||, quantifies the error in approx¬ 
imating the original tensor with the sum T. The 
number of components R is interpreted as the num¬ 
ber of sought mesostructures in the decomposition 
of the temporal network. Its value needs to be set 
by balancing the conflicting goals of minimizing e, 
i.e., recovering as much as possible of the original 
tensor, and avoiding overfittingl^^^IIlll sei;s 

of vectors C{i^ 2 ,...,r} can be 

combined into matrices A S B € 

and C G These matrices all have R columns, 

one for each component of the decomposition, i.e., 
one for each sought mesoscale structure. 

For an undirected temporal network (such as the 
one of our case study) the adjacency matrix repre¬ 
sented on each tensor slice is symmetric and the de¬ 
composition yields A ~ B. In practice it is possible 
to directly seek a decomposition with A = B, i.e., 
with component tensors of the form Sr = apOarOCr. 
Hence in the main text we only refer to the matri¬ 
ces A and C: the matrix elements air of A asso¬ 
ciate each component r = 1,..., i? with individual 
nodes of the network, while the matrix elements Ckr 
of C encode the temporal activity pattern of each 
component, with time indexed by k. The decom¬ 
position of Eq. can be written as 

R 

'Rijk — ^ ^ Rir ajr Cj^r • ( 6 ) 

r—1 

The matrix elements air indicate notion of member¬ 
ship of node i to component r, so this representa¬ 
tion is suitable to describe overlapping mesostruc¬ 
tures and complex correlations between connectiv¬ 
ity patterns and activity patterns over time. Notice 
that even when the original temporal network is 
unweighted, that is, the tensor T is binary-valued, 
the component tensors Sr correspond in general to 
weighted networks, and the approximating tensor 
T is in general real-valued. 


To compute the decomposition of Eq. we need 
to solve an optimization problem that minimizes 
the residual e. A standard way to do so is to con¬ 
vert the 3-mode problem yielded by Eq.j^into three 
coupled 2-mode sub-problems: This is done by un¬ 
folding the original tensor T along its three modes, 
a technique also known as matricizationl^. Using 
this representation, the original tensor approxima¬ 
tion problem is cast into 3 coupled matrix approx¬ 
imation problems that involve the factor matrices 
A, B, and C. Since the corresponding optimiza¬ 
tion problems are convex with respect to either of 
the matrices - but not with respect to them all - 
it is possible to solve the optimization through a 
method known as Alternative Non-negative Least 
Squares!^. The minimization of e is usually car¬ 
ried out with non-negativity and/ or sparsity con¬ 
straints on the factor matrice^^^llll^ as this is known 
to yield decompositions that can be interpreted as 
parts-based representations of the original data^^. 

Here we carry out the decomposition with 
non-negativity constraints, using Alternating Non¬ 
negative Least Squares and a block-coordina te de- 
scent method to improve convergence speed^^sMl 
Our implementation builds on the MATLAB code 
of ref. The number of components R is set 
by using the core consistency diagnosticPi^^ opti¬ 
mized over 5 stochastic realizations of the decom¬ 
position that differ for the random initial condi¬ 
tions of the factor matrices. The components Sr 
are ordered so that the Euclidean norm of Sr de¬ 
creases for increasing r. That is, the first column 
of matrices A, B and C correspond to the compo¬ 
nent with highest Euclidean norm, and successive 
columns encode weaker and weaker components. 
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Supplementary Information 


400 


Properties of the decomposition of the con¬ 
tact network 


The original temporal network 'T is decomposed by 
the non-negative tensor factorization method de¬ 
scribed in the Methods section and approximated 
by the resulting tensor T, which is a sum of R prod¬ 
ucts of lower-dimensional factors. This approxima¬ 
tion was carried out for several values of the number 
R of components, i.e. of rank-1 tensors, in order to 
determine the range of values such that R is large 
enough to describe the mesoscale properties of the 
original network and small enough to avoid overfit¬ 
ting. The core consistency, a measure introduced 
inl^ and plotted in Fig. SI for several values of the 
number R of rank-1 tensors, guides the choice of 
a representative number of components for T. In 



Figure SI: Core consistency curve with re¬ 
spect to the number of components. Each 
cross corresponds to a decomposition, 5 realizations 
are shown for each number of components used. 
The line acts as a guide for the eyes. 

the case study we consider, we limit the decompo¬ 
sition to i? = 14 terms to avoid overfitting. In Fig. 
|S2[ we show that this number of terms is sufficient 
by comparing the total activity (number of inter¬ 
actions in the whole data set) of each node in the 
approximated network R and in the original net¬ 
work T. As mentioned in the Methods section, T 
is real-valued even if T is binary-valued. 


Spreading processes and intervention strate¬ 
gies 

Both Susceptible-Infectible and Susceptible- 
Infectible-Recovered processes are run on the 
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Figure S2: Cumulative node strength in the 
original temporal network described by T, versus 
cumulative node strength in the reconstructed ten¬ 
sor T with i? = 14 components. 


reconstructed network and the altered networks. 
We consider i? = 14 altered networks: the r—th 
one, , is built by excluding the r—th component 
from the sum defining T. 

Details of the epidemics simulations 

The quantities measured for the SI process are aver¬ 
aged, for each altered network and for each parame¬ 
ter A, over all possible seeds and for each seed over 
50 spread starting times taken between 10% and 
30% of the total time length of the dataset. The 
quantities measured in the case of the SIR process 
are averaged, for each altered network and each set 
of parameters (A, /r), over all possible seeds and for 
each seed over 10 starting times taken between 10% 
and 15% of the total time length of the dataset. 


Efficiency of targeted interventions for the SI 
process: effect of different propagation proba¬ 
bilities 


Figure 2 of the main text shows the delay ratio of 
SI processes obtained by the selected removal of 
each mesostructure of the temporal network, for a 
specific value of the spreading rate A. The maxi¬ 
mal value of A tested is chosen such that A times 
the maximal weight of the links (here l.I) does not 
exceed I. 

We show in Fig. S3 the results obtained with dif¬ 
ferent values of A on each altered network together 
with the results obtained on 20 realizations of asso¬ 
ciated null models. Each null model associated to 
an altered network is built by taking the fully re- 
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constructed network T and zeroing out link weights 
at random until the sum removed is equal to the ac¬ 
tual cumulated weights removed in the correspond¬ 
ing altered network T’’. Although the efficiency of 
the single removal of each mesostructure depends 
on A, the global picture remains unchanged: the 
removal of one of the mixing patterns active during 
the first day (r = 12 or r = 14) largely outperforms 
the removal of any of the other mesostructures as 
well as the random removal of an equivalent amount 
of activity. 

Let us also note for precision that the half¬ 
infection time used in the computation of the infec¬ 
tion delay ratio is measured starting from the time 
at which the seed infects a node instead of the time 
of its first contact (as it was done in Ref. [M|) . The 
reason is that the links of the reconstructed net¬ 
work are not binary (active or inactive) but carry 
weights, so that the precise definition of a contact 
is blurred and recovering a binary definition of con¬ 
tacts would require to set an arbitrary threshold. 


Targeted intervention for the SIR process 

To quantify the difference between SIR processes 
unfolding on top of T and on top of altered net¬ 
works, we measured the averaged ratio, p\.^{r), be¬ 
tween the average size 11'' of the epidemic on the al¬ 
tered network and the average size 12 of the epi¬ 
demic on the full network T. To estimate the statis¬ 
tical significance of our results, we have performed 
Wilcoxon signed-rank tests for each parameter set 
(A, /i, r) comparing the epidemics size distributions 
pairwise between simulations performed on each al¬ 
tered network and on the full T. The test’s result 
shows that either the averages of the epidemic sizes 
are significantly different with a p—value lower than 
10“^, or the ratio p\,^{r) lies between 0.99 and 1. 
Therefore, this test shows that whenever p\^^{r) is 
smaller than 0.99, the observed difference between 
the outcomes of the spreading processes run on T'’ 
and T is significant. 

'We show in Fig. the ratio of epidemic sizes 
P\,fi{T') for SIR processes simulated on top of al¬ 
tered networks with a mesostructure r removed and 
on top of the whole decomposition T, for the values 
of r not shown in the main text. These mesostruc¬ 
tures describe activities occurring in single classes, 
and their removal have only negligible impact on 
the outcome of the SIR process. 


Evolution of the epidemic size in different 
regimes 


The temporal evolutions of the fraction of infected 
nodes are displayed in Fig. S5 for two sets of (A, p) 
chosen to be in two different regimes. In Fig. |S5[a), 
spread and recovery are slow, and the spread is not 
over at the end of the dataset. The removal of any 
of the mixing patterns thus has a similar impact 
at the end of the two days, independently of the 
fact that their activity occurs on the first or on the 
second day. For the parameter values used in Fig. 


S5 T) on the other hand, the spread is too fast to 


be mitigated. 
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Figure S3: Measured delay ratios for SI processes simulated on top of the reconstructed 
network T and of the 14 altered networks. The boxplots correspond to the delay ratios obtained 
by simulating SI processes on top of T and of 14 ensemble of networks altered by random removal of 
interactions: for each mesostructure r, we remove from T weights at random under the constraint that 
the sum of the removed weights is equal to the sum of the weights of the component r. Each boxplot 
gives the distribution of delay ratios obtained for 20 realizations of the random removal procedure. 
Horizontal lines inside each box show the median values, and each box extends from the first to the 
third quartile of the corresponding distribution. The reach of whiskers are determined by the sum of 
the third quartile and the interquartile weighted by 1.5. Points beyond the whiskers are outliers. The 
transmission probability of the SI processes is (a) A = O.I, (b) A = 0.3, (c) A = 0.5, (d) A = 0.9. 
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Figure S4: Ratio of epidemics sizes p\^^{r) for SIR processes simulated on top of altered 
networks with a mesostructure r removed and on top of the whole T. The values of r selected 
here are those which were not on Fig. 3 of the main text. The white part of each heat map corresponds 
to the region where the simulated epidemic process dies out. The black contour corresponds to the 
boundary of the parameter region such that the SIR process is over within the finite span of the dataset 
(i.e., to parameter values yielding epidemic processes that last at most two days): in this parameter 
region, the finite length of the dataset has no influence on the outcome of the epidemic. 
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Figure S5: Temporal evolution of the average fraction of nodes in the infected states for 

SIR processes simulated on top of altered networks with a mesostructure r removed and 

on top of the whole decomposition T. The values of r selected here are the same as in Fig. 4 of the 
main text, for (a) A = 2.10“^ and /i = 4.10“^ and (b) A = 1.6 10“^ and ^ = 4.10“^. The curve marked 
as none corresponds to the spread on the unmodified temporal network T. 
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